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SUMMARY 

This paper presents a solution of the problem of the 
theoretical flow of a frictionless incompressible fluid 
past airfoils of arbitrary forms. The velocity of the 2- 
dimensional floxo is explicitly expressed for any point 
at the surface, and for any orientation, by an exact 
expression containing a number of parameters which are 
functions of the form only and which may be evaluated by- 
convenient graphical methods. The method is particu- 
larly simple and convenient for bodies of streamline 
forms. The results have been applied to typical airfoils 
and compared vnth experimental data. 

INTRODUCTION 

The theory of airfoils is of vital importance in aero- 
nautics. It is true that the limit of perfection as 
regards efficiency has almost been reached. This 
attainment is a result of persistent and extensive 
testing by a large number of institutions rather than 
of the fact that the important design factors are known. 
Without the knowledge of the theory of the air flow 
around airfoils it is well-nigh impossible to judge or 
interpret the results of experimental work intelligently 
or to make other than random improvements at the 
expense of much useless testing. 

A science can develop on a purely experimental basis 
only for a certain time. Theory is a process of sys- 
tematic arrangement and simplification of known facts. 
As long as the facts are few and obvious no theory is 
necessary, but when they become many and less simple 
theory is needed. Although the experimenting itself 
may require little effort, it is, however, often exceed- 
ingly difficult to analyze the results of even simple 
' experiments. There exists, therefore, always a ten- 
dency to produce more test results than can be digested 
by theory or applied by industry. A large number of 
investigations are carried on with little regard for the 
theory and much testing of airfoils is done with insuffi- 
cient knowledge of the ultimate possibilities. This 
state of affairs is due largely to the very common belief 
that the theory of the actual airfoil necessarily would 
be approximate, clumsy, and awkward, and therefore 
useless for nearly all purposes. 

The various types of airfoils exhibit quite different 
properties, and it is one of the objects of aerodynamical 
science to detect and define in precise manner the fac- 



tors contributing to the perfection of the airfoil. 
Above all, we must work toward the end of obtaining 
a thorough understanding of the ideal case, which is 
the ultimate limit of performance. We may then 
attempt to specify and define the nature of the devia- 
tions from the ideal case. 

No method has been available for the determination 
of the potential flow around an arbitrary thick wing 
section. The exclusive object of the following report 
is to present a method by which the flow velocity at 
any point along the surface of a thick airfoil may be 
determined with any desired accuracy. The velocity 
of the potential flow around the thick airfoil has been 
expressed by an exact formula, no approximation hav- 
ing been made in the analysis. The evaluation for 
specific cases, however, requires a graphical determi- 
nation of some auxiliary parameters. Since the airfoil 
is perfectly arbitrary, it is, of course, obvious that 
graphical methods are to some extent unavoidable. 

Curiously enough, the theory of actual airfoils as 
presented in this report has been brought into a much 
simpler form than has hitherto been the case with the 
theory of thin airfoils. In the theory of thin airfoils 
certain approximations have restricted its application 
to small cambers only. This undesirable feature has 
been avoided, and the results obtained in this report 
have a complete applicability to airfoils of any camber 
and thickness. 

The author has pointed out in an earlier report that 
another difficulty exists in the theory of thin airfoils. 
It consists in the fact that in potential flow the velocity 
at the leading edge is infinite at all angles except one. 
This particular angle at which the theory actually 
applies has been defined as the ideal angle of attack. 
In the present work we shall not go any further into 
tins theory, since it is included in the following theory 
as a special case of rather limited practical importance. 

THEORY OF THICK AIRFOILS 

In the theory of functions there is a theorem by 
Eiemann 1 which shows that it is always possible to 
transform the potential field around any closed con- 
tour into the potential field around a circle. The 
direct transformation of an airfoil into a circle may, 

' Handbuch der Fhysik, Band III, p. 245, Fundanientalsatz der konformen 
Abbildung. 
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for analytical purposes, conveniently be performed in j 
two steps. The first step is to transform the airfoil into 
a curve which ordinarily does not differ greatly from a , 
circle by the transformation j 



, or 



(i) 



where f is a complex quantity defining the points in the | 
plane describing the flow around the airfoil and z ' \ 
is another complex quantity defining the points in the 
plane describing the flow around the almost circular ' 
curve. The constant a is of dimension length and is 
merely a geometrical scale factor. In the following 
theory, attention is directed to the fact that the shape 
of the curve resulting from transformation (I) is arbi- 
trary, since the airfoil shape is arbitrary. At a later 
point we shall transform this curve into a circle. 

The s' and the f planes are shown superposed in 
Figure 1. It will be noticed that at great distances 

z' 



This relation may further be conveniently expressed 
in hyperbolic functions 

j = 2a cosh <p cos 0 + 2ia sinh sin 8 

Since % = x + iy, the coordinates of the airfoil {x, y) are 
given by 

x = 2a cosh f cos d , TT > 
y = 2a sinh 4> sin 6 

We obtain a relation between 0 and the coordinates 
of the airfoil as follows: 

cosh 



2a cos 0 




sinh — -. — 7. 
r 2a sin 9 



and since cosh 2 ^ — sinh 2 $—1 
x 



( x Y-f v Y-i 

\2a cos dj \2a sin 0/ ' 



or developed 



where 



2sm 2 e = p + ^ 2 + gy 



(in) 



Similarly we obtain a relation between \j/ and the 
coordinates of the airfoil by using the equation 



Figure 1.— Showing the transformation from a noncircular curve B into an airfoil. 

from the origin z'—. that is, both flows are similar ! or developed 
at infinity. In particular, the "angle of attack," j 
defined as the direction of flow at infinitity with 
respect to some fixed reference line in the body, is 
identical in both flows. Near the origin the two flows 
are entirely different; one value of z' is, however, 
uniquely associated with a given value of f by the 
relation (I). 

We shall, at a later point, determine the flow in the 
z' plane. At present we shall determine the appear- 
ance of the airfoil when the almost circular curve B 
is given, or what amounts to the same tiling, we shall 
determine the curve B when the airfoil is given. In 
Figure 1, C is a circle of unit radius. Since the matter 
of dimensions is rather important, we shall avoid 
confusion in the following by adhering to this length 
as unit}-. The curve B is uniquely given by the rela- 
tion z' = ae^ +i6 where ^ is a known or unknown real 
function of the angle 0 where 6 varies from zero to 2it 
and i is the imaginary unit. Since the airfoil surface 
corresponds to the surface of the curve, the former is 
given from relation (I) as 



( * Yif y V 

\2a cosh \pj \2a sinh ty) 



Z=a(e* + e-*) cos d + i a(^-e~*) sin 0 



2 sinh 2 f = -p + tJp 2 + (|J (IV) 

Since f is generally small for wing sections it may be 
more conveniently expressed for purposes of calcu- 
lation as a series in terms oi^-. — ^> as follows: 

2a sin 0 

We have 

= sinh \j/ + cosh . 
= sinh f -f Vl + sinh 2 ^ , 

= 1 +sinh l/' + i sinh 2 if/+ ■ • • • 
<A=log» ^1 + sinh iA -t-^ sinh 2 f+ • • 
= sinh <P—^ sinh 3 \p+ • • • 

Kadb)'- 1 --"-- (IVa) 

[for,K log e 2] 



V 



2a sin 0 6 



We are now in a position ■ o reproduce the conformal 
representation of an airfoil in the r' plane, since tir 
each point of the airfoil (x, y) both 0 nnd ^ have been 
determined. 
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The curves ^ = constant are ellipses in the f plane 

7__^_Y+r y Y-i- 

\2a c ->sh £ / \2a sinh 
The foci are located at ( ± 2a, 0). The radius of eurva- 

, „ ; . . ■ (2a sinh 
ture at the end of the major axis is p = 2a cosh ^ 



or 



2a' 



(sinh i/Q 2 
• cosh^ 



(for small t//) 



Tliis relation is useful for the determination of f near 
the nose and the tail. 
The leading edge, corresponding to 0 = 0, is located at 



2a cosh 



<A 2 



= 2a + ai/' 2 = 2a-^p 



z 



Thus we see that the length 4a corresponds to the 
distance between the point midway between the nose 
and the center of curvature of the leading edge to the 
point midway between the tail and the center of 
curvature of the trailing edge. 2 

To establish the magnitude of the velocity at any 
point (:r, y) on the airfoil, we start in customary manner 
with the velocity around a circle in 2-dimensional 
flow. Contrary to usual practice we will, however, 
make the radius of the circle equal to ae^ where ^ 0 
. is a small constant quantity. This quantity is shown 
later in this report (equation (e)), to represent the 
average value of \p taken around the circle C. 

The potential function of the flow past this circle is 



(V) 



(VI) 



(reference 1, p. S3) and the velocity 3 

az \ , 2- / 2irz 

where- r is the circulation. This expression must 
vanish at the rear stagnation point 4 (Kutta condition) 
whose coordinate is z= — ae*< l+<ta+, r > i where a is the 
angle of attack and e T is shown to be the angle of 
zero lift. 

* The choice of axes is entirely arbitrary. It is a matter of convenience only to 
choose the axes so that the airfoil appears as nearly elliptical as possible, thereby 
making the "almost circular" curve B as nearly circular as possible by means of 
the single transformation I. It will be seen that the evaluation of the important 
integral appearing in the appendix is then most easily accomplished. In fact, the 
transformation I itself is only a matter of convenience to permit the ready evaluation 
of this integral. 



actually equals u—iv, the image of the velocity vector about the r-asi; 



' It is worthy of mention to note that the theory outlined in this report may actually 
be applied to smooth bodies of arbitrary shape if the circulation is specified. The 
term "wing sections" has been used in the title to imply bodies with sharp (or 
nearly sharp) trailing edges, whose circulation is or may be considered fixed by the 
Kutta condition or some equivalent assumption. 

149900—33 16 



We obtain T — **v(l-?£) 



e 



2i 



) 



=4irVae#° sin(a + € r ) (VII) 

This flow around the circle u,dj now be transformed 
into the flow around any other body. In the particular 
case in which the flow at infinity is not altered the 
circulation will not be altered and the force experienced 
by a body at the origin will remain at the fixed value 

i= P yr. 

We will now transform this circle, defined as 
z=ae fo+ip into our curve B denned by the relation 
3' = ae* +f9 . For this purpose we employ the general 

transformation z' = zen( A »+< B » ) j« which leaves the flow 
at infinity unaltered, the constants being determined 
by the boundary conditions. By definition 



Consequently 



or 



2 1 

2 1 

^ - $o + i (6 - <p) = n {A n + %B n ) ^ (cos n<p - i sin n <i) 

where z has been expressed in polar form 

2 = r(cos 'tp+i sin <p) 
and by De Moivre's theorem 

11, . . , 

-s = -s (cos n<p-i Bin n<p) 



Equating the real and imaginary parts we obtain the 
two Fourier expansions: 



and 



Bn 



'] 



(a) 
(b) 



e -v = nl7? cosn^-^- sin n^J 
A B 

The values of the coefficients ~> as well as the 
quantity fa, may be determined from (a) as follows: 



\[/ cos ntpd<p 



and 



r" irjo 
3 1 P 2t 

^ = -J 0 lMinn<pd P 

1 C 2r 
*«-2S Jo *** 



(c) 
(d) 

(e) 
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The quantity 6 — <p is necessary in the following 

A B 

analysis. Let us eliminate the coefficient ~~ and ~ 

in (b) by means of (c) and (d). 
We obtain 

2 1 f 2r 

(8—<p)c = ~, cos ?}><;- ipsmri(fd<9 
'* * Jo 

i r 2 ' 

-sin nif e - *pcosn<pd<p 
t Jo 

The subscript c is added to indicate that the angles so 
distinguished are kept constant while the integrations 
are performed. The expression may be simplified 

(6—<p) e =- n \ f (sin ncp cos n<p c — cos n<p sin n<p c )d</> 
But 

1 f , cos (2^+1)^=^ 



Hence 

dr 

d 



If , o 2 1 / , <r'\ 1 , , . _ , „. 
t' = 1 ~ = 2 "\ s ~ ?J = ? ^* a«-*"») 

= p [a(e*-c"*) cos «-fia (c*4-r#) sin 6} 
1 

= ■> [2a sinh $ cos 0-f 2{a cosh sin 5] 
Using the relations (II), 



2a sinh ^=-^— a and 2a cosh ^ = -^— a 
sm 6 T cos & 



we obtain 



d^* 1 

fig'**g/ (V cot e~ix tan 0). 



(IX) 



J 

It now remains to find the ratio jp> . From the 
relation 



we obtain 



sm 



<p — <Pc 



Therefore, 



or 



I cos (2/1+1) 
~2 Jo 



(v~ a c ) 



(*> _ <Pc) 



d<? But 



The latter integral is identically zero. (See Wilson, 
E. B. Advanced Calculus, p. 368. Follow method of 
exercise 10.) 
Then 



from which, 



dz' ,ri , d s, . , . _ . in 

= 2'^(^ + i + log Z) 

2 = aeV i,? 



2 = cLi (log 2) = d^ (log ° + + = Js ^ 



Therefore 



% =s 'ol^ + i (*-«?)+**>) 



For purposes of calculation this integral is expressed 
in convenient form in the appendix. 

We shall now resume the task of determining the 
velocity at any point of the surface of the airfoil. 

dw 



This expression may be written 



The velocity at the surface of the circle is 



.dz 



(see equation (VI) and footnote). For corresponding 
points on the curve B in the z' plane and on the airfoil 

in the £ plane the velocities are respectively • 



i But we have 



and 



dw dz dz' 



i or 



dz dz' df 

The quantities f and z' are related by the expression 



dz' , d , . , ... d$ 

-dz =s de^^'dz 

1 _ .dw 
z l dz 



--id <p=i d(& — d) + i dd 



and 
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where 



or 



da_' 
dz 



11.4 3 Cl J . , , 

dz' s' 



1 



dj 



dz 



1 + 6' 



(X) 



where e' and ^' indicate j| and -j~> 
Equations (IX) and (X) give now 



respectively, 



df 

3? 



d£' 

ds' 



^=-* (y cot 04iz tan 9) 



■ (y cot 04i.r tan 0) 



l 

3 14e' 



(XI) 



i = -VV cot 2 0 + ar tan- 0) (l+l /2 ) 
dz\ ae*° (l + O 



(XIa) 



The quantity c °t 2 0 + Qo) tan ^ * s rea< ^^ y 

to be equal to (by relation (II)) 

( o 0 Y+siu 2 9 
\2a sin 0/ 



or also 
Hence 



If 
i dz 



sinh 2 ^ + sin 2 9 



(1 4 e') 



(Xlb) 



jdw 

! dz , 



= F 2 [4 sin 2 (aieyJiS sin (a 4 e r ) sin Qx-r.<p) 
+ 4 sin 2 (a 4- <?)] 

' = 2 V[sin (a -f pi 4- sin (a -f e r )] 



Because we are interested more in the magnitude than, 
in the direction of the velocity we will write for the 
numerical value of this expression 



seen 



The numerical value of the velocity at the surface of 
' the circle is obtained by equations (VI) and (VII) as 
follows : 

Substituting the general point 3 = a^ 0+£(a+<?) , where 
a is the angle of attack as measured from the axis of 
coordinates, in equation (VI) 

^=-V(l-e-^«+^)-2i Fsin (o+ e T ) e -'^ 
= -V[L-C0s2(a4-p)42sm(a4e r )sm (* + <p) 
-ri (sin 2 (a4p)42 sin (a4e r ) cos (cc-r<p))] 



Eeplacing <p by 04 e <i T , the angle of zero lift, is the 
value of <p— 6 at the c \i? >> we have 

g^ =2 F[sin («404e)4sin (a + e r )] 
For a point on the airfoil we have, then, 

v = \ j • j~ ;' and from (XI), finally 

[Sin (a 4 0 + e) + sm (a±er)l£l_±_ e')^ 
V V A /(s"inhV + sin 2 e) (l+>' 2 ) {2UL) 

where tlie various symbols have the following signif- 
icance : 

v is the velocity at any point (a;, y) of the airfoil. 
Vis the uniform velocity of flow at infinity. 
y is the ordinate of the airfoil as measured from 
the z-axis, where to fix the system of coordi- 
nates (2a, 0) is the point midway between 
nose and center of curvature of the nose, 
and ( — 2a, 0) is the point midway between 
the tail and center of curvature of the tail. 
a is the angle of attack as measured from the 

x-axis as indicated in Figure 6. 
y, 0, yp, V, e, and e' are all functions of x. 
Equation (XII), expressing the value of the velocity at 
any point of an airfoil of any shape, is surprisingly simple 
when the complex nature of the problem is considered. 
It has the distinct advantage of being exact; no approxi- 
mations have been made in the preceding analysis. 

We shall note some of the properties of this impor- 
tant relation. Because y is generally small, the term 

=— a — -. is of influence chiefly near the leading edge, 
2a sm 0 J & & j 

where sin 0 is small. It is noticed, however, that if 

y = 

sin 0 

V= co. 

comes infinite for sinh ^=0 (thin airfoils). This fact 
has been pointed out in an earlier report. (Reference 

2.) The quantity 2 a ^m~6 01 ^ lS t * 1US °^ COn " 

siderable significance in the theory of thick airfoils. 

The velocity near the tail is obtained by putting 
0 = ir+A0 and e=e T + e'A6. Where A0 is a small angle, 
in equation (XII) 



= 0 for 0 = 0, equation (XII) yields in all cases 
This means that the velocity at the nose be- 



v\ 



we get 



= e t°Jl± !_' U sin (g + «+<0+sin (a + tr )] 
. v /(sirnrf4sin 2 0j (1 + 

e*> (1 + e') [-A04«4e r 4e'A04a4e r ] " 
" V(^ 2 + A0 2 ) (1 + *''*) 

e+» (14 V) 2 A0 
" A0 2 ) (1 4 iP 2 ) 

• (1 4V) 2 



(f) 



2 ) 



* It should be pointed out that the rear stagnation point is chostn to be on the 
i-aiis at 0=ir. The curvature at the tail is, as far as the specification of the ideal 
circulation is concerned, to be considered as a mechanical imperfection. 
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ip near the tail may be expressed as 



or 



The quantity is infinite if to is different from zero 

at A0 = O. The velocity is in this case zero, indicat- 
ing the presence of the rear stagnation point. If, on 
the other hand, to is zero, that is, if the tail is per- 
fectly-sharp, 

-f a =t' for A0 = O 
Ad 

and the velocitj r at the tail is 

+ e')' 

tT ~ 1 (i+toT 



or 



(g) 



(For the Clark Y, v T 2 is about 0.88 V 2 near the tail.) 
We obtain the front stagnation point by letting 
o = 0 in equation (XII). Hence 

a + 6 + e N = — (a + e T ) 

Q=-{2a + e N +e T ) 

In a previous report (reference 2) 

e jV + e T 
■ 2~ 

has been defined as the ideal angle of attack. It is 
seen that, for this angle of attack, 8 is zero or the stag- 
nation point occurs directly at the nose. 

Equation (XII) may also be applied to strut forms, 
and for such symmetrical shapes takes even a simpler 
form. 

PRACTICAL APPLICATION OF RESULTS 

We will now apply Formula (XII) to the typical 
case of the Clark Y airfoil and calculate the velocities 
at points of the airfoil surface. The detailed method 
of procedure is as follows. 

1. The axis of coordinates is drawn through the 
points (2a, 0) and ( — 2a, 0) located respectively at the 
point midway between the nose and the center of 
curvature of the nose and the point midway between 
the tail and the center of curvature of the tail. (See 
fig. 6.) The radius of curvature at the leading edge is 
1.75 per cent chord. 

2. The points {x, y) of the upper and lower surfaces 
of the airfoil are determined with respect to this axis. 



3. Sin 2 9, sin 6, and 6 are determined by the relation 
2 sin 'e-V + jjT^ -herep=l-(0-(fJ 

4. $ is given by the relation 

v \2a sin 6J Q\2a smdj 



6\2a 

5. \p is plotted as a function of 6 



1 P T 1 f 2 



6. Determine e c 



1 p* 
2t Jo 



f cot 



(<p—<Pc) 



dip by for- 



mula shown in the appendix: 

«„■=-- [0.62Sto c + 1 .065 (to - to-i) + 0.445 (to - to- 2 ) 
+ 0.23 1 (to - 4<-z) + 0. 104 (to - to-4)] 
where \p'c is the slope of the tp curve at <p — <p c to the 
value of yp at <p = <p c + -^, to at <p=<p e + ~jr , etc. 

to-i the value of \p at <p = <Pc~^, etc. 

7. From the e versus 5 curve and from the \p A'ersus 
6 curves e' and \p' are determined. 

8. Determine F by the relation 



(l + e')e* 



VlV 2 * sm 0 



y 



sin 



! 0~|(l + to 2 ) 



9. (0-f e) is determined in radians and degrees. 

10. Sin(0 + a+e)-f sin(a + e r ) is now calculated 
; where a is the angle of attack as measured from the 

axis of coordinates. 

11. y=F-[sin (<? + «+ e)+ sin {ct+e T )] 

■ P / » Y 

12. -=l-(j/j (pressure) 

i 

' The entire calculation, properly arranged, can be quite 
! accurately obtained in a very short time. 

COMPARISON WITH EXPERIMENTAL RESULTS 

In order to compare the theory with experimental 
I results, the geometric angle of attack a a as measured in 
' the wind tunnel must be corrected for a number of 
! items, such as finite span and effect of wall interference. 
! We may, however, obtain approximately the apparent 
j or effective angle of attack a A (in radians as measured 
j from the angle of zero lift) by taking the quotient 
; of the area of the pressure-distribution curve and 
] 5.5, since it is known that this value of the lift coeffi- 
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Figure 2.-Pressure-disfribation curves along x-aris of Clark Y; — against per cent chord 
(o) a=9°33'. (6) a=5° 19'. (c) a— 1" 10'. (<i) a=-3° 15' 
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FiGiritE 3.— Theoretical pressure distribution along y-a^is of Clark Y 
(a) a=9° 33'. (6) a = o° If. CO a=-l° IS'. (<f) o=-3° 15' 



cient is very nearly realized in most cases. This has 
been done in Table III, and the angle of attacks, which 
should be substituted in the Equation (XII), is 
given in the last column. The pressure distribution 
curves. Figures 2a, b, c, d, and 3a, b, c, d, were obtained 
by application of Equation (XII) to the Clark Y airfoil. 
Numerical results are shown in Tables I, II, and III. 
The experimental values are from original data sheets 
for N. A. C. A. Technical Report No. 353, and are 
not entirety consistent due to difficulties experienced 
in these experiments. After the theoretical pressure 
distribution curves have been obtained, the moments 
about any required axis may be found. Table IV 
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Figure 4.— Moment against angle of attack 

gives some of these results and Figure 4 shows the 
comparison with experimental data taken from N. A. 
C. A. Technical Report No. 312. 
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APPENDIX 



EVALUATION OF THE FORMULA 
= -± J%cot d*> 

Although the above integrand becomes positively 
and negatively infinite around <?=<? c , it. is readily 
verified that for $ finite, throughout the range 0 — 2tt, 
the integral remains finite, the positive and negative 
infinite strips exactly canceling each other. 

The value of the integral for any point <p c may be 
accurately obtained by the following device. We 
know that if f is a continuous function and the range 
tpi to <p 2 not too large • 

<P2~ <Pc 



i r 2 , t (<p-<p< 



sin 



dp is very nearly \p A log 



sin 



2 

where $ A is the average value of $ in the range <p x to 
<pi. Also near <p= <p t we may write 

W.+ fa-n) y t +( v - 9e yi^!+ .... 

Then for s a small quantity 

d*= = 2 f t 
= 4#' f 

(Since the even powers drop out and the lim <pcot a= 1). 

Let us now divide the interval 0 — 2tt into 10 parts, 
starting with <p c as a reference point. (See fig. 5.) 



pVot^ 

JVi-3 £ 



(f-ft). cot (e=30 dff 




Figure 5.— The ^against ^ curve, illustrating method of evaluation or u 



^ c ~10 Vc 




io t0 ^ 


, 3ir , 

^To'*- 


3w 

fo to 




JO tO ft~ 


7tt 

jTj'ft- 


7w 
jO to 


9tt 

io' 


Pc- Jq tO ft 


~i0'* e ~ 


9tt + 

io t0 ^ 




fo to 


oir ' 


jO to ft~ 


fjand 


3tT 

ft -10 to 





Then, 



1 f 2 * 



, («>-ft) , „ 
cot 2 ""V 



1 



sin: 



3tt 



J + log 

20 



sin 



. 57T . 77T 

sin 20 an go 

t (fa - is) log 3- + (h ~ is) log =- 



+ (.ii~i-i) log 



. 3ff 
sin 20 

sin 20 

7-7T 



sm 



20 



sm 



20J 



= -- [0.62S i' c + 1.065 + 0.445 (i 2 ~is) 

IT 

+ 0.231 + 0.104 (&-^_ 4 )] 

where is the slope of the ^ curve at <p= ft 



11 value of i at <p — <p c + = » at u>= <s c — ±> 

o 0 

1 2 at ¥>= ft-!--^' ^3 at p = ¥><!+ -J- » e tc. 

To evaluate the above integral it is, strictly speak- 
ing, necessary to know ^ as a function of <p rather than 
of 0. 1 We have tp = 6 + e. For all flattened or stream- 
line bodies, however, e is small; for ordinary airfoils 
it is, in fact, so small that \p(6) may unconditionally be 
considered equal to i(4>). For the sake of mathe- 
matical accuracy we will, however, indicate how the 
problem may be solved also for bodies of more irregular 
contour by successive approximations. We have 

*G»W W + •*'(*) + • • • 

As a first approximation we neglect the second and all 
following terms of this expression. The value of « 
. thus obtained by graphical integration or otherwise 
is- then used in the expression for i(<p) an d a second 
integration is performed, etc. 



1 The equation for e is a nonlinear integral equation and to obtain its exact solu- 
tion is a difficult matter; fortunately because of the small magnitude of £ the solu- 
tion is obtainable to any desired accuracy by ordinary definite integrals. 
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APPLICATION OF FLOW FORMULA TO THE SPECIAL I 
CASE OF AN ELLIPTIC CYLINDER 

As a matter of interest we will assume the form of j 
the body to be the ellipse (&^)' 1 | 

and find y for zero angle of attack, i. e., we have ' 

yp = io = constant, i£' = 0, e=0, e' = 0, a = 0. 
Equation (XII) becomes 

sin Q-e* 



■(f)" 



y-e* 



2a sinh f 



and 



Vsinh 2 l/' + sin 2 6 I . , , , , / 

V Smh 'K2asinh*. 



y 



P. 
i 



■ay 



,2* 



(2a2/) 2 



(2a smb. </-)*+ (2 cn/) 2 



This result checks exactly with the form given by 
Dr. A. F. Zahm in N. A. C. A. Technical Report No. 
253, Flow and Drag Formulas for Simple Quadrics, 
equation 14. 
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EXPLANATION OF THE TABLES 

The first part of Table I refers to the upper surface 
or to positive ordinates of the Clark Y, the second 
part to the lower surface or to negative ordinates. 
Column 1 gives the location in per cent of the chord; 2 
gives the ordinates with respect to the z-axis in this 
same unit; 3 and 4 give x and y in the present sj r stem 
of coordinates; 5, 6, and 7 give sin 2 d, sin 0, and 0, 
respectively (Equation (III)) ; S gives $ (by equation 

(IYa)); 9 gives e (appendix); 10 and 11 give -J^ and ^ 

as obtained from \p against 6 and e against 6 curves; 
(See figs. 7 and 8). Column 12 gives the quantity 



(1 + e') e*» 



(e*'=i:il); 



Column 14 gives 6 + e in degrees. The velocity at any 

point x and angle of attack a is given b} r 

i-=F[sin (a-r0-!-e)+ sin (a + e T )]-F and the pressure, 

It must be noted that a is measured from the line of 




Figure 6.— Clark Y airfoil— showing system of coordinates 

flow to the z-axis as shown in Figure 6, and if otherwise 
measured, must be reduced to this basis. 




Tail 



Figuee 7.— The unit circle 2=e'*, the circle z=e* rfi *, and the corresponding 



curve z'—e" 




Upper si .-fjee 

e 

Fiouse 8— (a) T„e & against 6 curve for the Clark Y. (6) The t against B curve for the Clark ^ 
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TABLE I 
CLARK Y 



UPPER SURFACE 





V in % • 


X 




sin'0 

- 


sin B 


8 ■ 
radians 


1> 


c 




i* 


F 


S+e 
radians 




0 


0 


2.035 


0.000 


0.000 


0.000 




0.000 


0.188 


-0.079U ff ) 


0.000 


0.075 


6.33 


-0. 079 


e 

4 


/ \ 
32 


1.25 


1.S9 


1. 985 


.0804 


.0474 


.218 


.220 • 


.184 


-.067 


.035 


.080 


4.20 


.153 


8 


47 


2.50 


3.08 


. 1.934 


.124 


.0994 


.315 


.321 


.193 


—.055 


.075 


.105 


3.28 


.266 


15 


17 


.5.0 


4.59 


1.833 


.185 


.195 


.442 


.458 


.208 


-.038 


.060 


.115 


2.52 


.420 


24 


4 


7.5 


5.61 


1. 732 


.226 


.282 


.531 


.559 


.213 


-.026 


.030 


.115 


2.16 


.533 


30 


33 


10 


6.45 


1.631 


.260 


.365 


.605 


.649 


.214 


-.016 


.020 


.120 


1.92 


.633 


36 


16 


15 


7.70 


1. 430 


.311 


.512 


.716 


.798 


.216 


.003 


.000 


.120 


1. 66 


.801 


45 


53 


20 


8.55 


1.228 


.345 


.640 


.801 


.929 


.214 


.020 


-. 030 


.120 


1.48 


.949 


54 


22 


30 


9.23 


.824 


.372 


.837 


.915 


1.156 


.202 


.044 


-.060 


.105 


1. 30 


1.200 


68 


45 


40 


9.28 


.421 


.374 


.957 


.978 


1. 361 


.190 


.063 


-. 065 


.091 


1.21 


1.424 


81 


35 


50 


8.74 


.0175 


.353 


1.000 


1.000 


1.571 


.176 


.080 


-.0S5 


;078 


1. 17 


1.651 


94 


35 


60 


7.77 


-. 386 


.313 


.963 


.982 


1.761 


.159 


.094 


-. 100 


.070 


1.18 


1. 855 


106 


16 


70 


6.32 


-.790 


.255 


.847 


.920 


1.972 


.137 


.106 


—.104 


.025 


1. 22 


2.078 


119 


3 


80 


4.48 


-1. 193 


.181 


.649 


.806 


2.204 


.112 


.110 


-. 107 


.004 


1.36 


2.314 


132 


34 


90 


2.40 


-1. 596 


.097 


.368 


.606 


2 490 


.080 


.106 


-.110 


-. 015 


1.78 


2.596 


148 


43 ■ 


95 


1. 25 


-1.79S 


.050 


.195 


.443 


2.6S3 


.056 


.095 


-.080 


—.075 


2.30 


2.778 


159 


10 


100 


.06 


-2.001 


.002 


.000 


.000 


3. 142 


.030 


.062( ir ) 


-.053 


-.080 


Large 


3.204 


183 


33 



LOWER SURFACE 



0 


0 


2.035 


0.000 


0.000 


0. 


ooo 


6.2S3 


0. 1SS 


-0. 079 


0.000 


0. 075 


6. 33 


6.204 


-4 


32 1 


1. 25 


1.53 


1. 935 


-.0617 


.0387 




197 


6. 085 


. 1S6 


-. 100 


.170 


.085 


4.71 


5.985 


-17 




2.50 


1. 95 


1.934 


-. 0787 


.0822 




287 


5. 992 


.137 


-. 105 


.162 


.055 


3.64 


5.887 


-22 


« 


5.0 


2.38 


1.833 


-.0960 


.171 




414 


5.857 


.116 


-.110 


.130 


.000 


2.55 


5.747 


-30 


43 I 


7.5 


2.61 


1. 732 


-.105 


.258 




508 


5. 750 


. 103 


-.107 


.110 


-.015 


2.09 


5.643 


-36 


40 1 


10 


2.73 


1.631 


-.110 


.342 




585 


5. 60S 


.094 


-.105 


.093 


-.032 


1.81 


5.553 


-41 


50 1 


15 


2.84 


1.430 


-. 115 


.492 




702 


5.505 


.0S2 


-.100 


.080 


-.040 


1. 49 


5.399 


-50 


40 


20 


2.78 


1.228 


-.112 


.623 


— 


791 


5.371 


.071 


-.093 


.075 


-. 04S 


1. 33 


5.278 


-57 


35 


30 


2.47 


.824 


-. 0996 


.831 




912 


5.135 


.055 


-. 081 


.058 


-.050 


1. 15 


5. 054 


-71 


0 


40 


2.12 


.421 


— . 0SS5 


.956 




973 


4. 922 


.045 


-.072 


.043 


-.050 


1.075 


4.850 


-82 


6 


50 


1.78 


.0175 


-. 0720 


1.000 


-1 


000 


4.712 


.036 


-.058 


.040 


-.060 


l.Ol 


4.154 


-93 


21 


60 


1. 38 


-. 3S6 


-. 0556 


.962 




981 


4.518 


.028 


-.015 


.030 


-.060 


1.06 


4.473 


-103 


40 


70 


1.03 


-.790 


-.0416 


.844 




919 


4.306 


.023 


-.034 


.025 


-. 060 


1.13 


4.272 


-115 


14 


SO 


.74 


-1. 193 


-. 0295 


.645 




803 


4.075 


.013 


-.0 f 


.022 


-.070 


1.275 


4.055 


-127 


39 


90 


.40 


-1.596 


-.0161 


.364 




604 


3. 790 


.013 


.010 


.018 


-.0S2 


1.685 


3. 790 


-142 


51 


95 


.24 


-1. 798 


- 0097 


. 191 




438 


3. 595 


.011 


.023 


-.017 


-.090 


2.30 


3.618 


-152 


43 


100 


.06 


-2.001 


-.002 


.000 




000 


3. 142 


.030 


.062 


-.053 


-.080 


Large 


3.204 


-176 


27 



TABLE II 
CLARK Y 
|fi=[sta (a-r-a+O+sin (a+«p] • F f =l"jff 
ct=9° 33' 

« T =3° 33' sin (a-i-.p =sin 13°6'=0.2267 



Upper surface 



Lower surface 



%c 


e+a 




! sin 
1 («+<*+<) 


sin 1 
(94-CZ-M + 
sin (a+« r ) 


1*1 

.V| 


I r ' 

v: 


P 
Q 




sin 
(9+b+O 


sin 
sin (ct+tj.) 


|»l 

1 v! 




P 

8 


0 


0 

5 


/ 
1 


1 

! 0. 0875 


0.3142 


1.99 


1 

3.96 


-2.96 


0 

5 


t 
1 


0. 0875 


0. 3142 


1.99 


3.96 


-2.96 


1.25 


18 


20 


■ .3145 


.5-112 


2.26 


5.12 


-4. 12 




32 


-. 1311 


.0956 


.450 


.203 


.797 


2.50 


24 


50 


.4200 


. 6467 . 


2. 12 


4.49 


-3. 39 


--li 


8 


-. 2272 


.0005 


.002 


.000 


1.000 


5.0 


33 


37 


! .5537 


.7801 


1.97 


3.87 


-2. 87 


-21 


10 


-.3611 


-. 1344 


.343 


.117 


.8S3 


7.6 


40 


6 


.6441 


.8703 


1.88 


3.55 


-2.55 


-27 


7 


-. 4558 


-. 2291 


.479 


.229 


.771 


10 


45 


49 


.7171 


. 9433 


1. 82 


3.30 


-2.30 


-32 


17 


-. 5341 


-. 3074 


.556 


.303 


.692 


15 


55 


26 


.8235 


1.0502 ■ 


1.74 


3.04 


-2.04 


-41 


7 


-. 6576 


-. 4307 


.643 


.414 


.586 


20 


63 


55 


. 8981 


1. 1248 


1.68 


2.83 


-1.83 


-48 


2 


-. 7435 


-. 5168 


.685 


.469 


.531 


30 


78 


18 


.9792 


1. 2059 


1.57 


2.47 


-1.47 


-61 


27 


-.8784 


-.6517 


.748 


.558 


.442 


40 


91 


8 


. 9998 


1. 2265 


1.49 


2 °2 


I 22 


-72 


33 


-.9540 


— . 7273 


.788 


.621 


.379 


50 


104 


8 


. 9697 


1. 1964 


1.40 


1.98 


-o!is 


-83 


43 


-.9941 


-. 7674 


.800 


.610 


.360 


60 


115 


49 


.9000 


1. 1267 


1.33 


1.76 


-0.76 


-94 


7 


-. 9974 


-. 7707 


.821 


.674 


.326 


70 


128 


36 


.7815 


1.0082 


1.23 


1.51 


-0.51 


-105 


41 


-. 9627 


-. 7357 


.834 


.696 


.304 


80 


142 


7 


. .6141 


. S40S 


1.14 


1.31 


-0.31 


-118 


6 


-. SS21 


—.6554 


.830 


.698 


.302 


90 


158 


16 


.3703 


. 5970 


1.06 


1.13 


-0. 13 


-133 


18 


-. 7278 


-.5011 


.844 


.711 


.289 


95 


168 


43 


. 1957 


.4224 


.97 


.93 


0.05 


—143 


10 


-. 5995 


-. 372S 


.857 


.735 


.265 


100 








! 










Formula (g) 




.935 


.875 


.125 



Table II gives the numerical values for Figure 2a in detail as an example. See also Table I. 



TABLE III 



TABLE IV 



Geometric 
angle as 
measured 

from 
chord line 
experi- 
mentally 



13 25 

7 25 

-1 40 

-4 35 



I 



Area of j 
eiperi- ; 
mental ' 
pressure : 
distri- 



5.5' 



bution ] radians 
curve ! 



1 
t 

Apparent 

angie 
measured 

from 
zero lift 

degrees 



1.260 
.855 
.221 
.030 



0.229 
. 155 
.040 
.0053 



13 6 

8 52 

2 17 

0 18 



Angle 
measured 

from 
chord line 
(to be used 
as a basis 
for 

comparison 
with a g ) 
i=a.-3°33' 



9 33 
5 19 



-3 15 



Figure 


Mm 

Moment 
about line 
z=25 
per cent 
chord 


M, 

Moment 
about line 
y=0 


Mr 

Moment 
about 
point 
x=25 

per cent 
chord, 
V=0 


2a 


-0. 0896 


-0.00S5 


-0.093 


2b 


-.098 


-.005 


-.103 


2c 


-.098 


.000 


-.098 


2d 


-.093 


.0015 


-.091 



Positive directions of axes and angles (forces and moments) are shown by arrows 



Axis 


Force 
(parallel 
to axis) 
symbol 

^ 


Moment about axis 


Angle 


Velocities 


Designation 


Sym- 
bol 


Designation 


Sym- 
bol 


Positive 
direction 


Designa- 
tion 


Sym- 
bol 


Linear 
(compo- 
nent along 
axis) 


Angular 


Longitudinal... 
Lateral 


X 
Y 
Z 


X 
Y 
Z 




L 
M 
N 


Y >Z 


roll.. - 


$ 
0 


u 

V 

■ . w 


P 
1 

■ r 


pitching 


Z — JT 




Normal 


X > Y 











Absolute coefficients of moment 
L „ M 



qbS 



qcS 



C '~gbS 



Angle of set of control surface (relative to neu^ 
teal position), 5. (Indicate surface by proper 
subscript.) - ; - 



4. PROPELLER SYMBOLS 



2> 



— :_, . r: " 



Diameter. • *- / \ r~ -2/ . ■ \ -{ 
Geometric pitch. 4 . 
p/D, Pitch ratio. 
V\ Inflow velocity. " 
V„ Slipstream velocity. 

T, Thrust, absolute coefficient C r =— ^ 
Q, Torque, absolute coefficient Cq= J^ps 

5. NUMERICAL RELATIONS 



P, Power, absolute coefficient C P = yptyf : " 

(7 S , Speed power coefficient =-yp-2« 

ii, Efficiency. • ^ 4 ~. - 
», Revolutions per second, r. p. s. 

*, Effective Tielix angle = tan" 1 



1 hp =76.04 kg/in/s=550 lb./ft./sec, 
1 kg/m/s = 0.01315 hp 
1 mi./hr. = 0.44704 m/a; .'" 
1 m/s=>2.23693 mi./hr." . . 



,1 lb. =0.4535924277^;^ 
1; kg =2,2046224 lb, 

lm= 3:2808333 ft. ' 



